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ABSTRACT 

We report the results of Chandra and XMM-Newton observations of a new outburst of an ultraluminous super- 
soft X-ray source in M101. CXOUJ140332. 3+542103 was observed in a low luminosity state (L x ~ 10 37 ergss" 1 ) 
between 2004 January and May. The low-state X-ray spectra were relatively hard; the combined low state spec- 
trum can be fitted with a combination of a power-law with photon index of 1 .4 and a blackbody of 63 eV. During 
2004 July, the source underwent a strong outburst and the peak 0.3-7 keV luminosity reached 3 x 10 40 ergs s , 
with a bolometric luminosity of about 10 41 ergs s -1 . The outburst spectra were very soft and can generally be fit- 
ted with a blackbody model with temperatures of 50-100 eV. In two of the observations, absorption edges at 0.33 
keV, 0.56 keV, 0.66 keV, and 0.88 keV were found. An XMM-Newton observation was also performed during the 
decay of the outburst, and a power-law tail was seen in addition to the supersoft spectrum. We consider different 
accretion models; one involving an intermediate-mass black hole can explain the observations. 

Subject headings: black hole physics — galaxies: individual (M101) — X-rays: binaries — X-rays: galaxies 



1. INTRODUCTION 

Recent high resolution X-ray observations reveal that there 
are large number of ultraluminous X-ray sources (ULXs) in 
many nearby galaxies. ULXs are usually defined as very bright 
(Lx ~ io 39-41 ergss" 1 ) off-nucleus, variable X-ray sources in 
galaxies. Since the luminosity exceeds the Eddington luminos- 
ity of a 10M Q black hole (BH), ULXs are considered to be can- 
didates for intermediate-mass black holes (IMBHs). More re- 
cently, a number of ULXs have been found to be supersoft; their 
X-ray spectra can be fitted with a blackbody model with tem- 
perature of less than 100 eV (e.g., Swartz et al. 2002; Mukai et 
al. 2003; Fabbiano et al. 2003; Di Stefano & Kong 2003; Kong 
&DiStefano 2003). 

In this Letter, we report a series of Chandra observations 
of the ultraluminous SSS in M101 (CXOU J140332.3+542103) 
from the low state (Kong & Di Stefano 2004) to outburst (Kong 
2004a, 2004b). Shortly after the discovery of the outburst, we 
initiated a multiwavelength campaign to monitor the source. 
We here also present the results of an XMM-Newton follow- 
up observation conducted about two weeks after the first out- 
burst observation. The primary result of radio observations by 
the Very Large Array was reported by Rupen, Sjouwerman, & 
Kong (2004). We will present the detailed analysis of the opti- 
cal and radio observations in a follow-up paper. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra 

M101 was observed by Chandra in 2004 January, March, 
May, and July. Table 1 summarizes the Chandra observations. 
We also list the two Chandra observations taken in 2000 March 
and 2000 October, and the follow-up XMM-Newton observation 
(see § 2.2). In all observations, the center of M101 and the ultra- 
luminous SSS are in the S3 chip of the ACIS-S. All data were 
taken in very faint mode (VFAINT). We reprocessed the raw 

1 http://asc. harvard, edu/ciao/threads/aciscleanvf/ 



data to utilize the VFAINT mode and to remove the 0.5-pixel 
randomization. However, for data taken during the outburst, we 
found that VFAINT processing flagged some source events of 
the SSS as background and we therefore did not apply VFAINT 
processing for the data taken in July 1 . In order to reduce the 
instrumental background, we screened the data to allow only 
photon energies in the range of 0.3-7 keV. We also searched 
for periods of high background using source free regions in the 
SI chip and filtered all the high background periods. Only the 
2004 January 24 and 2004 May 7 observations were affected; 
the exposure times listed in Table 1 correspond to the good time 
intervals. All data were reduced and analyzed with the CI AO 
v3.1 package, and calibration database CALDB 2.27. 

2.2. XMM-Newton 

Shortly after the discovery of the outburst of 
CXOUJ140332.3+542103 with Chandra, we requested XMM- 
Newton director's discretionary time for a follow-up observa- 
tion. The XMM-Newton observation was taken on 2004 July 23 
for 30 ks, and the instrument modes were full-frame, medium 
filter for the three European Photon Imaging Cameras (EPIC). 
After rejecting those intervals with a high background level, we 
considered a good time interval of 23 ks. Only data in 0.2-12 
keV were used for analysis. Data were reduced and analyzed 
with the XMM-Newton SAS package v5.4.1. 

3. ANALYSIS AND RESULTS 

We performed spectral analysis for all Chandra data taken 
between 2004 January and 2004 July. The two datasets taken 
in 2000 were reported in Mukai et al. (2003) and Di Stefano & 
Kong (2003). Spectral analysis was performed by using Sherpa. 
We also used XSPEC vl 1 .3 for independent check. In each ob- 
servation, we extracted the source spectrum from a 6" radius 
circular region centered on the source, while an annulus region 
centered on the source was used as the background. During the 
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outburst, there are enough photons to rebin the spectra with at 
least 15 counts per spectral bin, and used \~ statistics to find the 
best-fitting parameters. For the data during the low state (Jan- 
uary, March, and May), very few source photons were received 
(~ 10-40) in each observation, we therefore performed spec- 
tral analysis with unbinned data (background not subtracted) 
for which C-statistics (Cash 1979) were employed. To examine 
the low state spectrum in detail, we also combined all the low 
state data to perform spectral fit with binned data and \ 2 statis- 
tics. We tried power-law, blackbody, and disk blackbody mod- 
els with interstellar absorption. Table 2 and 3 show the best- 
fitting spectral parameters for all Chandra observations taken 
in 2004, and the X-ray spectra are shown in Figure 1 . 

Table 1 
Observation Log of M101 



For the low state data, power-law models provide better fits 
than blackbody models. Each low-state observation is con- 
sistent with an absorbed power-law model; the photon index 
ranges from ~ 0.8-2.6 among the observations, with an aver- 
age of 1 . 8 . The Nh converged to zero and we therefore fixed it at 
the Galactic value (1 .2 x 10 20 cm" 2 ; Dickey & Lockman 1990). 
The 0.3-7 keV luminosity is about (1.5-4.6) x 10 37 ergss _1 . 
Fixing the Nh at a higher value (1.5 x 10 21 cm" 2 ; see below), 
the luminosity is about a factor of 2 higher. The combined low 
state data, however, cannot be fitted with any single-component 
model. A simple power-law model produces a soft excess. 
We therefore fitted the spectrum with a composite blackbody 
and power-law model (Nh = 1.5 x 10 21 cm -2 , kT = 63 eV, and 
a = 1.4; see Table 2); the power-law component contributes 
roughly 70% (30%) of the absorbed (unabsorbed) 0.3-7 keV 
flux. 

During the 2004 outburst, blackbody models provide good 
fits to the spectra. Disk blackbody model is also equally ac- 
ceptable; the inner disk temperatures are slightly higher than 
the blackbody temperatures, but the normalizations (hence the 
inner disk radii) often have large error ranges. In general, the 
blackbody temperatures are between ~ 50- 100 eV, and the 
Nh is about (1-4) x 10 21 cm" 2 , significantly higher than the 
Galactic value. We also fitted the spectra with the Nh fixed at 
1.2 x 10 20 cm" 2 . However, all the fits are unacceptable with 
reduced \ 2 between 1.5 and 2.4. The 0.3-7 keV luminosity 
ranges from 10 39 ergss _1 to 10 40 ergss _1 . We also derived the 
bolometric luminosity from the normalization of the blackbody 
model and it is between 5 x 10 39 ergss _1 and 10 41 ergss" 1 . In 
two cases, the spectra are complex and single component mod- 
els cannot provide good fits to the data. We included two ab- 
sorption edges in each case and the fits were significantly im- 
proved. Specifically, we found absorption edges at 0.66 keV 
and 0.87 keV on 2004 July 6 but the edges were not observed 
on 2004 July 8. Instead, we found two new edges at 0.33 keV 
and 0.56 keV. 

We used similar techniques to analyze the XMM-Newton 
data. We extracted the source spectrum from a 15" radius 



circle centered on the source. The source was fainter than 
previous Chandra observations with about 160 background- 
subtracted counts in EPIC-pn. Since the source was detected 
only in EPIC-pn, we fitted the EPIC-pn spectrum with several 
single-component models. However, all the fits are unaccept- 
able. We therefore fitted the spectrum with a composite black- 
body and power-law model. The model consists of a supersoft 
component (53 eV) and a hard power-law tail (a = 0.7); the 
power-law component contributes about 40% (8%) of the 0.3- 
7 keV absorbed (unabsorbed) flux. The 0.3-7 keV luminosity 
is 6. 1 x 10 38 ergs s" 1 , indicating that the source was in the decay 
stage of the outburst. 
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FIG. 1. — Unfolded spectra of the ultraluminous SSS. The models corre- 
spond to the fits listed in Table 2 and 3. The total spectrum, blackbody com- 
ponent, and power-law component are shown in green, red, and blue, respec- 
tively. 

4. DISCUSSION 

Although over 100 ULXs have been discovered (Miller & 
Colbert 2003), only a handful of sources have L x > 10 41 ergs 
s , assuming isotropic radiation (Matsumoto et al. 2001; Gao 
et al. 2003; Davis & Mushotzky 2004; Soria & Motch 2004). 
A luminosity this high is difficult to achieve in an X-ray bi- 
nary unless the accretor has a mass greater than roughly 10M Q . 
We do have evidence that CXOU J140332.3+542103 is likely 
to be an X-ray binary, since its luminosity has been observed 
to change by a factor of > 10 on a time scale of hours (Mukai 
et al. 2003; Di Stefano & Kong 2003). While it is possible that 
our luminosity estimates are higher than the true luminosity, the 
effects that lead to overestimates, such as beaming, or various 
anisotropies, tend to change estimated luminosities downward 
by a factor of roughly 10. Its unusually high bolometric lumi- 
nosity (1.2 x 10 41 ergs s" 1 ), coupled with its short-time-scale 
time variability, therefore make CXOU J140332. 3+542103 a 
good candidate for an accreting IMBH. For instance, King et 
al. (2001) proposed that anisotropic X-ray emission can result 
super-Eddington luminosity for a stellar-mass BH. However, 
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Table 2 

Best-fitting Spectral Parameters for the Low State 
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NOTES — All quoted uncertainties are 90% confidence. Sources without reduced \ 2 an d dof values were fitted with unbinned data using CASH statistics; 
"goodness-of-fit" determined by Monte-Carlo simulations (MC) is employed and parameters of both blackbody and power-law models are shown. 
" in units of 10 21 cm~ 2 . 

* For low state data, we list the probability that the best fit model would produce a lower value of the CASH statistics than that calculated from the data, as 

determined via XSPEC Monte-Carlo simulations. A best fit model should have a value of about 0.5. 

c 0.3-7 keV luminosity (x 10 38 ergs s" 1 ), assuming 6.7 Mpc (Freedman et al. 2001). 

d 0.3-7 keV luminosity (x 10 38 ergs s" 1 ) derived from the combined power-law and blackbody model. 

e fixed 

f Combined low-state data modeled by an absorbed power-law and blackbody model. 



in order to achieve such a high luminosity for a stellar-mass 
BH, extreme beaming is required and the disk is expected to be 
much hotter. Similarly, the temperature of radiation pressure- 
dominated accretion disk model proposed by Begelman (2002) 
is too high. The pure blackbody spectrum also makes relativis- 
tically beaming (Kording et al. 2002) unlikely. 

Its soft spectrum in the high state is another important piece 
of evidence. An accretion disk around a very massive BH is ex- 
pected to produce supersoft X-ray emission (see Di Stefano & 
Kong 2002). For instance, if the supersoft emission comes from 
the inner disk, a 70 eV SSS at luminosity of 10 40 ergss" 1 cor- 
responds to a ~ 6000 Mq BH. If we used the 90% lower limits 
of the inner disk temperature derived from the disk blackbody 
fits, the BH mass is estimated to be > 2800M Q (Makishima et 
al. 2000). This is consistent with prior work on IMBH models 
for ULXs (see e.g., Miller & Colbert 2003). It also comple- 
ments other work on evidence for cool disks in ULXs that has 
also been considered as evidence for IMBH models (Miller et 
al. 2003,2004) although a power-law component usually con- 
tributes a significant fraction of X-ray emission. Furthermore, 
the high state luminosity is approximately 16% of the Edding- 
ton luminosity for a 6000 M Q BH; we therefore expect the inner 
disk to be optically thick, which is consistent with the IMBH 
interpretation. 

The state changes of CXOU J140332.3+542103 span 
a luminosity range larger than a factor of 1000. 
CXOU J 140332.3+542 103 was first detected by ROSAT in 1996 
November (Wang et al. 1999), with an extrapolated unabsorbed 
0.3-7 keV luminosity of 8 x 10 38 ergss" 1 (blackbody model 
with N H = 10 21 cm -2 , kT = 75 eV). Higher absorption (3 x 10 21 
cm" 2 ) would imply a luminosity up of 4 x 10 39 ergss" 1 . The 
source was not detected in other /?(95Arobservations (see Wang 
et al. 1999). In a 108 ks ROSATHRl observation taken in 1996 
May/June, we can set a 3a detection limit at 6.8 x 10 37 ergs s" 1 
(power-law model with Nh = 1.2 x 10 20 cm" 2 and a = 2). Set- 
ting the Nh at 1.5 x 10 21 cm" 2 , the luminosity limit becomes 
1.5 x 10 38 ergss" 1 . In 2000 March and October, Chandra 



detected the SSS in the very soft state with luminosities of 
- 10 39 ergss" 1 (Mukai et al. 2003; Di Stefano & Kong 2003). 
In 2002 June, XMM-Newton observed M101 and did not de- 
tect the SSS (Jenkins et al. 2004). We re-analysed the image 
and found that the SSS may be barely detected with possi- 
ble contamination from nearby sources. We derived the 0.3-7 
keV luminosity of the SSS as ~ 10 37 ergss" 1 , by assuming a 
power-law model. The source was also in the low state between 
2004 January and May with luminosity of ~ 10 37 ergss" 1 . It 
is clear that the SSS has had at least three major outbursts 
(L x > 10 39 ergss" 1 ). On the other hand, the very low lumi- 
nosities (~ 10 37 ergss" 1 ) during 2002 and 2004 indicate that 
the source was in the low state. The source varies by as much 
as a factor of 1000 between the low state and the high state. 
This amplitude is even greater than many Galactic BHs. Re- 
markably, the source also shows spectral changes. It is very 
clear that there is a power-law component (a = 1.4) in the 
composite low state spectrum, while the high state spectra are 
supersoft. Our XMM-Newton observation taken during the de- 
cay of the outburst revealed that the blackbody component was 
still strong and there was a very hard power-law tail, similar to 
the low state. While it is difficult to make direct comparison 
with Galactic BHs because of the blackbody emission in the 
low state and the pure blackbody spectrum in the high state 
(see however Zdziarski & Gierlinski 2004), it is unusual for a 
SSS to have high/very soft state and low/hard state transition. 
With a peak bolometric luminosity w 10 41 ergs s" 1 , it is likely 
that IMBH is the central engine of the system. Interestingly, 
the blackbody component is always seen while the power-law 
component becomes stronger when the source is at lower lu- 
minosities. It may indicate that the power-law component is 
due to Comptonization of soft photons (Wang et al. 2004). The 
photon index is, however, harder than that of Galactic BHs in 
the low state (a ~ 1.7). None of the Galactic BHs has similar 
spectrum (McClintock & Remillard 2003). The closest exam- 
ple is Galactic microquasar V4641 Sgr for which the photon 
index was measured to be between 0.6 and 1.3 (Miller et al. 
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Table 3 

Best-fitting Blackbody Parameters for the High State 
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2002; McClintock & Remillard 2003). 

Another important feature of our X-ray spectral fits is the 
presence of absorption edges, which we found at 0.33, 0.57, 
0.66, and 0.88 keV in two of the high state spectra. These 
are consistent with C V, N VI, N VII, and O VIII edges. We 
note, however, that the 0.33 keV edge may be due to calibration 
of the ACIS-S near the carbon edge around 0.28 keV. These 
features may signal the presence of highly ionized gas in the 
vicinity of the accretor (e.g., warm absorber), and may be con- 
sistent with an outflow from the source. In fact, outflow mod- 
els have been suggested for ULXs. King and Pounds (2003) 
suggested that the soft component of some ULXs is due to 
outflow of a stellar-mass BH. Mukai et al. (2003) also used 
similar argument to explain the super-Eddington luminosity of 
CXOUJ140332.3+542103. However, the new data which re- 
veal the changes of temperature and bolometric luminosity, the 
extremely high luminosity, and the state transition of the source 
will be difficult to explain by such a model. It remains a puzzle 
that we saw different edges in the two observations, possibly 
related to the geometry of the system. Nevertheless, outflows 
from IMBHs may be expected. 

It is also interesting to note that similar absorption edges are 
expected and have been observed in white dwarf (WDs) sys- 
tems such as the recurrent nova U Sco (Kahabka et al. 1999), 
CAL 87 and RX J0925.7^1758 (Ebisawa et al. 2001). SSS 
radiation with luminosities in the range from 10 37 ergs s" 1 to 
roughly the Eddington limit for a 1.4M Q WD are expected 
for nuclear burning WDs (van den Heuvel et al. 1992). In- 
deed, such models have been proposed to explain SSSs in 
the Galaxy and Magellanic Clouds. Note, however, that neu- 



tron star and BH models have also been proposed. The high 
state luminosity of this source is almost two orders of mag- 
nitude larger than that of any other known SSS, and it rules 
out steady nuclear burning on hot WD models as an explana- 
tion for CXOU J140332.3+542103. A similar argument can be 
made for neutron star models. A nova explosion could explain 
some features of the data. It is not favored, however, because 
the photospheric radius and lack of hard radiation at the peak 
of the outburst have no obvious explanation. 

It may be impossible to conclusively establish, using cur- 
rent technology, that any X-ray source in an external galaxy is 
an IMBH. A system like CXOU J140332.3+542103 is there- 
fore particularly valuable, because it provides 4 different 
pieces of evidence that together make a consistent argument 
in favor of an IMBH interpretation. Its high-state luminos- 
ity, its short-time-scale variability, its soft high-state spec- 
trum, and its pronounced spectral changes, all suggest that 
CXOU J140332.3+542103 is a strong IMBH candidate. 
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